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Abstract 

Up to now attempts to measure the general relativistic Lense- 
Thirring effect in the gravitational field of Earth have been performed 
by analyzing a suitable J2 — J4— free combination of the nodes f2 of 
LAGEOS and LAGEOS II and the perigee oj of LAGEOS II with the 
Satellite Laser Ranging technique. The claimed total accuracy is of 
the order of 20-30%, but, according to some scientists, it could be an 
optimistic estimate. The main sources of systematic errors are the 
mismodelling in the even zonal harmonic coefficients J; of the multi- 
polar expansion of the gravitational potential of Earth and the non- 
gravitational perturbations which plague especially the perigee of LA- 
GEOS II and whose impact on the proposed measurement is difficult 
to be reliably assessed. Here we present some evaluations of the accu- 
racy which could be reached with a different Ji— free observable built 
up with the nodes of LAGEOS and LAGEOS II in view of the new 
preliminary 2nd-generation Earth gravity models from the GRACE 
mission. According to the GRACE-only based EIGEN-GRACE02S 
solution, a 1-sigma upper bound of 4% for the systematic error due to 
the even zonal harmonics can be obtained. In the near future it could 
be possible to perform a reliable measurement of the Lense-Thirring 
effect by means of the existing LAGEOS satellites with an accuracy of 
a few percent by adopting a time span of a few years. The choice of a 
not too long observational temporal interval would be helpful in reduc- 
ing the impact of the secular variations of the uncancelled even zonal 
harmonics J4 and Jg whose impact is difficult to be reliably evaluated. 
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1 Introduction 



The general relativistic gravitomagnetic force induced by the gravita- 
tional field of a central rotating body of mass M and proper angular mo- 
mentum J is still awaiting for a direct, unquestionable measurement. Up to 
now there exist some indirect evidences of its existence as predicted by the 
General Theory of Relativity (GTR in the following) in an astrophysical, 
strong-field context [2j and, in the weak-field and slow-motion approxima- 
tion valid throughout the Solar System, in the fitting of the ranging data to 
the orbit of Moon with the Lunar Laser Ranging (LLR) technique 3 . The 
measurement of the gravitomagnetic Schiff precession of the spins of four 
spaceborne gyroscopes jlj in the gravitational field of Earth is the goal of 
the Stanford GP-B mission [3] which has been launched on April 2004. The 
obtainable accuracy should be of the order of 1% or better. 

The Lense-Thirring effect on the geodesic path of a test particle freely 
falling in the gravitational field of a central rotating body [Hj consists of 
tiny secular precessions of the longitude of the ascending node f2 and the 
argument of pericentre u of the orbit of the test particle 
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where a, e and i are the semimajor axis, the eccentricity and the inclination, 
respectively, of the orbit, c is the speed of light and G is the Newtonian 
gravitational constant. 

The LAGEOS III/LARES mission [3|H] was specifically designed in order 
to measure such effect, but, up to now, in spite of its scientific validity 
and relatively low cost, it has not yet been approved by any space agency 
or scientific institution. Recently, a drag-free version of this project [S], 
in the context of the relativistic OPTIS mission is currently under 
examination by the German Space Agency (DLR). 

2 The current LAGEOS-LAGEOS II Lense-Thirring 
experiment 

Up to now, attempts to observationally check the Lense-Thirring effect in the 
gravitational field of Earth have been performed by analyzing the accurately 
recovered orbits of the existing LAGEOS and LAGEOS II satellites with the 
Satellite Laser Ranging technique (SLR) ^2^3- The adopted observable 
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is 

5Q L + 0.295<«7 L 11 - 0.35<fo L 11 ~ 60.2^ LT , (2) 

where 50 and 6u are the orbital residuals of the rates of the node and 
the perigee and /^lt is the solved-for least square parameter which is in 
Newtonian mechanics and 1 in GTR. The Lense-Thirring signature, entirely 
adsorbed in the residuals of ft and to because the gravitomagnetic force has 
been purposely set equal to zero in the force models, is a linear trend with 
a slope of 60.2 milliarcseconds per year (mas yr _1 in the following). The 
standard, statistical error is evaluated as 2%. The claimed total accuracy, 
including various sources of systematic errors, is of the order of ^2 El 
20 — 30%. However, such estimate would be too optimistic according to 
some scientists who propose a different error budget |13j . 

2.1 The error budget 

The main sources of systematic errors in this experiment are the unavoidable 
aliasing effect due to the mismodelling in the classical secular precessions 
induced on Q and cu by the even zonal coefficients J/ of the multipolar expan- 
sion of geopotential and the non-gravitational perturbations, which sever ly 
affect the perigee of LAGEOS II O [El HE1 El HH1 , whose impact on the 
proposed measurement is difficult to be reliably assessed. It turns out that 
the mismodelled classical precessions due to the first two even zonal har- 
monics of geopotential J2 and J4 are the most insidious source of error for 
the Lense-Thirring measurement with LAGEOS and LAGEOS II. The com- 
bination © is, by construction, insensitive just to J2 and J4. According to 
the full covariance matrix of the EGM96 gravity model the error due 
to the remaining uncancelled even zonal harmonics amounts to almost 13% 
|20| (1-sigma calculation). However, if the correlations among the even zonal 
harmonic coefficients are neglected and the variance matrix is used in a 1- 
sigma Root-Sum-Square fashion, the error due to the even zonal harmonics 
of geopotential amounts to 46.6% [201 - Such approach is considered more 
realistic by some authors ^3] because nothing assures that the correlations 
among the even zonal harmonics of the covariance matrix of the EGM96 
model, which has been obtained during a multidecadal time span, would 
be the same during an arbitrary past or future time span of a few years 
as that used in the LAGEOS-LAGEOS II Lense-Thirring experiment. A 
1-sigma upper bound of almost 83% for the gravitational error can be ob- 
tained by adding the absolute values of the individual errors 20 . Another 
point to be emphasized is that the use of the perigee of LAGEOS II forces 
to adopt an observational time span of many years in order to view certain 
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long-period harmonic perturbations of gravitational and non-gravitational 
origin as empirically fitted quantities which can be removed from the time 
series without corrupting the extraction of the genuine relativistic secular 
trend. Indeed, it turns out that the perigee of LAGEOS II is affected by the 
ocean tidal perturbation Ki, 1 = 3, p = 1, which has a period of 5 years 
|21| . and by a direct solar radiation pressure harmonic with a period of 11.6 
years |14j . According to |14| I15j . the non-gravitational part of the error 
budget would amount to 28-30% over seven years. Moreover, recent reex- 
aminations of certain nonconservative accelerations acting upon LAGEOS 
II would suggest that it could be reduced down to a 13% level over the 
same time span. However, it must be pointed out that such estimates are 
based on certain refinements of the non-gravitational force models which 
were not included ^Ej m the GEODYN II orbit processor used at the time 
of the analysis of especially as far as certain tiny non-gravitational per- 
turbations of thermal origin j£] are concerned. Moreover, it must also be 
recognized that the estimates of the authors of ^3] are different from such 
evaluations; indeed, it can be argued that their evaluation of the impact of 
the nonconservative accelerations on the measurement of the Lense-Thirring 
effect with the perigee of LAGEOS II reported in is of the order of 48- 
99%, if the optimistic 13% error, based on the EGM96 full covariance, is 
adopted. 

3 The role of the new Earth gravity models from 
the CHAMP and GRACE missions 

From the previous considerations it could be argued that, in order to have 
a rather precise and reliable estimate of the total systematic error in the 
measurement of the Lense-Thirring effect with the LAGEOS satellites it 
would be better to reduce the impact of geopotential in the error budget 
and/or discard the perigee of LAGEOS II which is very difficult to handle 
and is a relevant source of uncertainty due to its great sensitivity to many 
non-gravitational perturbations. 

The forthcoming more accurate Earth gravity models from CHAMP |23| 
and, especially, GRACE [21] will yield an opportunity to realize both these 
goals, at least to a certain extent. 
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3.1 The EIGEN-GRACE02S model 



In order to evaluate quantitatively the opportunities offered by the new 
terrestrial gravity models we have preliminarily used the recently released 
EIGEN-GRACE02S gravity model It is important to note that such 
model represents a long-term averaged GRACE-only solution (110 days); 
moreover, the released sigmas of the spherical harmonic coefficients of the 
geopotential are not the mere formal statistical errors, but are calibrated, 
although preliminarily. Then, guesses of the impact of the systematic error 
due to the geopotential on the measurement of the Lense-Thirring effect 
based on this solution should be rather realistic. However, caution is advised 
in considering the so obtained evaluations because of the uncertainties of the 
calibration process which affect especially the even zonal coefficients |13j . 

With regard to the three-elements combination ([%]). it turns out that 
the systematic error due to the even zonal harmonics of the geopotential, 
according to the variance matrix of EIGEN-GRACE02S up to degree I = 70, 
amounts to 0.2 mas yr _1 (Root Sum Square calculation), yielding a 1-sigma 
0.4% error in the Lense-Thirring effect. The sum of the absolute values of 
the individual errors yields an error of 0.4 mas yr _1 , i.e. a 1-sigma upper 
bound of 0.7% in the Lense-Thirring effect. Of course, even if the LAGEOS 
and LAGEOS II data had been reprocessed with the EIGEN-GRACE02S 
model, the problems posed by the correct evaluation of the impact of the 
non-gravitational perturbations on the perigee of LAGEOS II would still 
persist. 

A different approach could be followed by taking the drastic decision 
of canceling out only the first even zonal harmonic of the geopotential by 
discarding at all the perigee of LAGEOS II. The hope is that the resulting 
gravitational error is reasonably small so to get a net gain in the error budget 
thanks to the fact that the nodes of LAGEOS and LAGEOS II exhibit a 
very good behavior with respect to the non-gravitational perturbations. 
Indeed, they are far less sensitive to them than the perigee of LAGEOS 
II. Moreover, they can be easily and accurately measured, so that also the 
formal, statistical error should be reduced. A possible observable is [22j 

5Q L + 0M65Q L 11 ~ 48.2/iLT- (3) 

A similar proposal can be found in [21], although numerical details are not 
released there. According to the variance matrix of EIGEN-GRACE02S up 
to degree / = 70, the residual signal due to the even zonal harmonics from 
I = 4 to / = 70 is 1.5 mas yr _1 (Root Sum Square calculation), i.e. a 
1-sigma 3% systematic bias in the Lense-Thirring effect. The sum of the 
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absolute values of the individual errors yields an upper bound of 1.9 mas yr 
, i.e. a 1-sigma 4% systematic error. EGM96 would not allow to adopt 
j3J) because its full covariance matrix up to degree I = 70 yields an error 
of 47.8% while the error according to its diagonal part only amounts even 
to 104% (1-sigma Root Sum Square calculation), with an upper bound of 
177% (1-sigma sum of the absolute values of the individual errors). Note 
also that the combination @ preserves one of the most important features 
of the combination of (J2J) of orbital residuals: indeed, it allows to cancel 
out the very insidious 18.6-year tidal perturbation which is a I = 2, m = 
constituent with a period of 18.6 years due to the Moon's node and nominal 
amplitudes of the order of 10 3 mas on the nodes of LAGEOS and LAGEOS II 
|21| . On the other hand, the impact of the non-gravitational perturbations 
on the combination (j3J) over a time span of, say, 7 years could be quantified 
in 0.1 mas yr -1 , yielding a 0.3% percent error. The results of Table 2 and 
Table 3 in |H] have been used. It is also important to notice that, thanks to 
the fact that the periods of many gravitational and non-gravitational time- 
dependent perturbations acting on the nodes of the LAGEOS satellites are 
rather short, a reanalysis of the LAGEOS and LAGEOS II data over just 
a few years could be performed. As already pointed out, this is not so 
for the combination (J2J) because some of the gravitational an d non- 
gravitational |14j perturbations affecting the perigee of LAGEOS II have 
periods of many years. Then, with a little time-consuming reanalysis of 
the nodes only of the existing LAGEOS and LAGEOS II satellites with 
the EIGEN-GRACE02S data it would at once be possible to obtain a more 
accurate and reliable measurement of the Lense-Thirring effect, avoiding the 
problem of the uncertainties related to the use of the perigee of LAGEOS 
II. 

The choice of an observational time span of just a few years would also 
be helpful in reducing the impact of the secular variations of the uncancelled 
even zonal harmonics J4 and Jq. The problem of the impact of the secular 
variations of the even zonal harmonics on the proposed measurements of the 
Lense-Thirring effect has never been addressed, up to now, in a satisfactorily 
way. In [23| it has been claimed, perhaps too superficially, that the secular 
variations of the zonals would not affect the combination © because they 
can be accounted for by an effective J| ff . In fact, the effective J2 means 
a lumped effect that it has not been possible to separate with one or two 
satellites. But the individual effects are still there; they just getted blurred. 
One can use the lumped effect to get some insight into the total error in 
the secular rates of the zonals, but it tells nothing about the individual 
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contributions. So, it is not possible to cancel them out in the combination 
(J3J). For the topic of J/, which has recently received great attention by the 
geodesists' community in view of unexpected variations of 2 J2, see [2Ej. By 
assuming 5 j 4 = 0.6 x 1CT 11 yr _1 and 5Jq = 0.5 x 10~ n yr _1 [25] . it turns 
out that the 1-sigma error on the combination Q would amount to 2.1% 
over one year. However, it must be pointed out that it is very difficult to 
have reliable evaluations of the secular variations of the higher degree even 
zonal harmonics of geopotential also because very long time series from the 
various existing SLR targets are required. 

3.2 Alternative combinations 

In |27| the possibility of using a multisatellite linear combination including 
the the nodes of LAGEOS, LAGEOS II, Ajisai, Starlette and Stella has 
been investigated. It is, by construction, insensitive to the first four even 
zonal harmonics of the geopotential. On the other hand, the inclusion of 
the nodes of the other existing SLR satellites, which orbit at much lower 
altitudes than the LAGEOS satellites, introduces, in principle, much more 
noise from the higher degree even zonal harmonics which such combination 
would be sensitive to. The hope was that the improvements in the knowledge 
of just the higher degree even zonal harmonics from the new GRACE-based 
solutions would make such an alternative combination competitive with the 
two- nodes combination of (jHJ). The recent results from EIGEN-GRACE02S 
rule neatly out this possibility. Indeed, it turns out that, if, on the one 
hand, the problem of the J; would be greatly reduced, on the other hand, 
the 1-sigma Root Sum Square percent error would be 13% with an upper 
bound of 42% from the sum of the absolute values. 

More favorable and interesting, at least in principle, is the situation with 
another node-only combination proposed in [22] . It includes the nodes of LA- 
GEOS, LAGEOS II, Ajisai and the altimeter satellite Jason-1 whose orbital 
parameters are similar to those of Ajisai; apart from the LAGEOS satellites, 
Ajisai and Jason-1 have the most interesting orbital configurations, among 
those of the existing accurately tracked satellites, for our purposes. The 
weighing coefficients of the nodes are 1 for LAGEOS, 0.347 for LAGEOS II, 
-0.005 for Ajisai and 0.068 for Jason-1; the gravitomagnetic slope is 49.5 mas 
yr _1 . It turns out that the EIGEN-GRACE02S model yields a systematic 
gravitational error of 1% (1-sigma Root Sum Square calculation) and an 
upper bound of 2%. Moreover, J4 and Jq would not affect this combination. 

X I am grateful to J. Ries for helpful discussions on this problem. 

2 Fortunately, any issues concerning J2 do not affect the combination 
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However, the possibility of effectively getting long time series of the node 
of Jason-1 should be demonstrated in reality. Finally, dealing suitably with 
the non-gravitational perturbations acting on it in a genuine dynamic way 
would be a very demanding task. 

4 Conclusions 

When more robust and complete terrestrial gravity models from CHAMP 
and GRACE will be available in the near future, the two-nodes/LAGEOS- 
LAGEOS II combination © could allow for a measurement of the Lense- 
Thirring effect with a total systematic error, mainly due to geopotential, 
of the order of a few percent over a time span of some years without the 
uncertainties related to the evaluation of the impact of the non-gravitational 
perturbations acting upon the perigee of LAGEOS II. The choice of a not too 
long observational time span should also be helpful in keeping the systematic 
error due to the secular variations of the even zonal harmonics below the 
10% level. 

On the other hand, the obtainable accuracy with the node-node-perigee 
combination (J2J), whose error due to geopotential will remain smaller than 
that of ©, is strongly related to improvements in the evaluation of the 
non-gravitational part of the error budget and to the use of time spans 
of many years. However, it neither seems plausible that the error due to 
the non-conservative forces will fall to the 1% level nor that a reliable and 
undisputable assessment of it will be easily obtained. 

Alternative combinations including the orbital data from the other ex- 
isting SLR satellites are not competitive with the combination of 

A combination including also the nodes of the SLR Ajisai and altimeter 
Jason-1 satellites, together with the nodes of the two LAGEOS satellites, 
would be slightly better from the point of view of the reduction of the 
systematic error due to the geopotential, in particular with regard to the 
effects of the secular variations of the even zonal harmonics. However, this 
gain could be lost due to the difficulties of dealing with the non-gravitational 
perturbations affecting the nodes of Jason-1. 
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